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ABSTRACT 



We present the results of high resolution mapping of the CO(1-0) and CO(2-l) emission of the ultraluminous infrared galaxy (ULIRG) 
IRAS 10190+1322, with the IRAM interferometer, down to an angular resolution of ~0.3". This object is composed of two interacting 
galaxies with a projected nuclear separation of 6kpc, and was selected to analyze the physical and dynamical properties of the 
molecular gas in each galaxy in order to study the conditions that lead a galaxy pair to become ultraluminous in the infrared. With the 
exception of Arp 220, the closest ULIRG, this is the first time that the CO emission is morphologically and kinematically resolved 
in the two interacting galaxies of a ULIRG system. In one of the galaxies the molecular gas is highly concentrated, distributed in 
a circumnuclear disk of 1.7 kpc in size. The molecular gas in the presumably less infrared luminous galaxy is distributed in a more 
extended disk of 7.4 kpc. The molecular gas mass accounts for ~10% of the dynamical mass in each galaxy. Both objects are rich 
enough in molecular gas, ~ 4 10 9 Mq, as to experience an infrared ultraluminous phase. 
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1. Introduction 

Ultraluminous infrared galaxies (ULIRGs: > 10 12 Lq) 
present signs of current or past interactions in their optical 
and near-infrared images in a fraction that approaches 100% 
dMurphv et al.lll996tlBushouse et al.ll2002tlVeilleux et al.ll2002|) . 
The ULIRG phenomenon has been observed along the full in- 
teraction phase, from widely separated pairs of galaxies (up to 
a projected nuclear separation Sn ~ 50 kpc) to advanced and 
relaxed mergers. However, the distribution of nuclear separa- 
tions is highly peaked at low values wit h 60% of the IRAS 
Uy ULIRG sample having S N < 2.5 kpc dVeilleux et al.ll2002l) . 
In particular, this has hindered the study of the molecular 
gas properties in ULIRGs, sub-arcsecond resolution observa- 
tions being necessary to morphologically and kinematically re- 
solve the emission of the two interacting galaxies. In fact, this 
has been achieved only in Arp 220, t he mos t nearby ULIRG 
dDownes & Solomon! Il998l hereafter |PS98t ISakamoto etafl 
119991) . 

To date, there have been two millimeter interferometric sur- 
veys devoted to study the molecular gas content of ULIRGs with 
separated pairs of galaxies, but they have lacked the necessary 
sensitivity and/or angular resolut ion to detect and resolv e the 
CO emission of the two galaxies. lDinh- V-Trung et al. (200l]) de- 
tected CO(1-0) emission with the BIMA array in five out of 
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Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
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For the definition of L„ see Table 1 in ISanders & Mirabell £l996). 
In this paper we will assume a flat A-do minated cosmology described 
by H = 71 kms- 1 Mpc~' and n m = 0.27 (Spergel et al. 20Q). 



six targets (20 kpc < 5n < 51 kpc), but only in one compo- 
nent of the ULIRG pairs. Interferometric observations with the 
O VRO array of five double nuclei ULIRGs (3 kpc < S n < 5 kpc) 
bv lEvans et ail d2002l) have succeeded in detecting the CO(1-0) 
emission of both nuclei in two cases, but did not resolve then- 
individual molecular gas distributions. The overall conclusion of 
all these observations, mostly derived from upper limits, was that 
>2/3 of the molecular gas mass lies around one of the nuclei, be- 
ing this nucleus the most active of the two (LINER or Seyfert op- 
tical classification). Additional observations are required to con- 
firm or reject this conclusion. The recent upgrade of the IRAM 
Plateau de Bure Interferometer (PdBI) to larger baselines cur- 
rently allows this sort of studies. 

IRAS 10190+1322 (L ir = 10 I2 L o at D L = 340 Mpc) 
is composed of two interacting galaxies with a projected nu- 
clear separation of 6 kpc (4"). The south-west (SW) galaxy 
is brighter than the north-east (NE) galax y in the optical but 
the latter is brighter in the near-infra r ed dMurphv et all [ 200 1 ; 
iKim et al] 120021; lim anishi etalj 120061: iDasvra et al] 120061) . At 
15 GHz iNagar et al] d2003l) detected compact radio continuum 
emission only in the NE galaxy. Th e system was not spa tially 
resolved at 4.86 GH z and 1 .4 GHz bv lCrawford et al] d 19961) and 
Whit e et all d 1997b . but the radio continuum distributions peak 
towards the NE galaxy position. All these observations point to 
the NE galaxy to host the bulk of the mid- and far-infrared lu- 
minosity of IRAS 10190+1322. However, limanishi et al.ld2007l) 
found that the Spitzer 5.2-14.5/im spectra of the NE galaxy is 
only s lightly brighter than that of the SW galaxy. iMurphv et ail 
(2001) obtained Paa images of both galaxies and concluded that 
the emission in the NE galaxy is more compact and absorbed 
than in the SW galaxy. This result agrees with their kinemat- 
ical analysis where the two galaxies are approaching to expe- 
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Fig.l. a) CO(1-0) integrated intensity map of IRAS 10190+1322 at the reference position {a nm , (5 J2 ooo) = (10 h 21 m 42 s :6, 
13°06'54'.'4). Contours for the CO integrated intensity are lcr, 2<r, 4<x, 6<x and 9<x to 39<x in steps of 6<x (lcr = 0.21 Jy kms" 1 
beam 1 ). The synthesized beam size, shown in the bottom left corner, is 1.22" x 0.57" with a position angle of 27°. The black 
crosses mark the central position of the NE and SW galaxies (see Table Q]) determined by elliptical Gaussian fits in the uv plane, b) 
CO(1-0) velocity map relative to the redshift of observation (z D b s = 0.07656). Contours for the CO velocity are -300 to 200km s _1 
in steps of 50 km s _1 . The straight lines indicate the position angle of the mayor axis in each galaxy estimated from the velocity 
channel maps, c) CO(2-l) integrated intensity map. Contours for the CO integrated intensity are lcr, 2<x, Act, 6cr and 9<x to 33<x in 
steps of 6cr (lcr = 0.45 Jy kms -1 beam -1 ). The synthesized beam size is 0.61" x 0.28" with a position angle of 21°. d) CO(2-l) 
velocity map. Contours for the CO velocity are -300 to 200km s _1 in steps of 50km s . In order to show the more extended and 
weaker emission coming from the SW galaxy, the moment maps have been derived applying a 3cr-clipping to each channel of the 
CO(1-0) and CO(2-l) clean maps in the velocity range between -360 and 240km s _1 . 



rience a second encounter after the first close passage, and the 
NE galaxy has resulted more affected by the interaction than 
the SW galaxy due to the initial configuration of the encounter. 
Rup ke et all (120021) classified the NE galaxy as a LINER and the 
SW galaxy as an HII galaxy. iFarrah etall d2003l) modelled the 
overall spectral energy distribution (SED) of the system and de- 
termined that the contribution to L u from an active galactic nu- 
cleus (AGN) could be as high as the 20%. llmanishi et all (120061 
120071) did not find evidences of an obscured AGN in the L-band 
and 5.2-14.5/mi spectra of both galaxies. 



2. Observations 

IRAS 10190+1322 (z = 0.07656) was observed with the six 
antennae of the PdBI in February 2006 in the new extended 
A configuration (baselines up to 760 m). The 3 mm and 1mm 
receivers were tuned to the CO(1-0) and CO(2-l) redshifted 
lines at 107.074GHz and 214.143GHz, respectively. The sys- 
tem temperatures were 140 K for CO(1-0) and 250 K for CO(2- 
1). Four correlator units covered a total bandwidth of 580 MHz 
at each frequency. Phase and amplitude calibrations were done 
on nearby quasars. 

Reduction using the GILDAS software provided data cubes 
with a resolution of 0.12"/pixel and 0.06"/pixel that were 
cleaned with the standard method with a velocity resolution 
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Fig. 2. a) IRAM 30 meter telescope CO(l-O) spectrum of IRAS 
10190+1322 complex (light grey). The CO(1-0) line profile of 
the NE galaxy observed with the PdBI is represented in dark 
grey at the same panel, b) The result of the substraction of the 
NE galaxy PdBI CO(1-0) line to the total single dish CO(l- 

0) line is shown in light grey. In dark grey we plot the PdBI 
CO(1-0) line spectrum of the SW galaxy, c) and d) PdBI CO(2- 

1) line profiles of the NE and SW galaxies, respectively. The 
PdBI spectra of each galaxy were derived integrating inside the 
area that includes all pixels with an emission >3cr in any of the 
30kms _1 channel maps. 

of 30kms~'. The restored clean beams are 1.22" x 0.57" (PA 
= 27°) at 2.8 mm and 0.61" x 0.28" (PA = 21°) at 1.4mm. 
The rms noise levels in the cleaned maps (at 30kms _1 veloc- 
ity resolution) are 1.57mJybeam _1 for the CO(1-0) line and 
4.25 mJy beam -1 for the CO(2-l) line. We tentatively detected 
2.8 mm continuum emission of 0.9 mJy, just below the 3<x level, 
towards the central position of the NE galaxy. This density flux 
is compatible with the 60 fim and 100//m IRAS density fluxes 
of the whole system assuming optically thin grey-body emission 
with a dust temperature T^ust ~ 38 K and a spectral index of the 
dust absorption coefficient ft ~ 1 .4. 

IRAS 1 1 90+ 1 322 was also observed with the IRAM 30 me- 
ter telescope in June 2006. We used 2 SIS 3 mm receivers to ob- 
serve the 2 polarizations simultaneously at 107.074GHz. At this 
frequency, the telescope half-power beam width was 23" and 
the antenna temperature to flux conversion factor was S/T* = 
6.3Jy/K. During the observations the typical system tempera- 
ture was 170 K (on the T* scale). The observations were done in 
wobbler switching mode, with reference positions offset by 2' in 
azimuth. Two 1 MHz filter banks provided a total bandwidth of 
512MHz, or about 1400 kms -1 . 

3. Results 

Emission of the first two CO rotational transitions has been de- 
tected and spatially resolved towards both galaxies. A summary 
of the observational results is given in Table Q] In Fig. Q] we 
have plotted the integrated intensity and velocity maps of the 
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Table l.IRAS 10190+1322 main properties 





Unit 


NE Galaxy 


SW Galaxy 


"JZUUU 




10 h 21 m 42i754(l) 


10 h 21 m 42!490(4) 






13°06'55'.'59(1) 


13°06'53'.'80(5) 


V — C 7«k* 
» *- <.obs 


km s 


-100(10) 


40(10) 


*CO( 1—0) 


Jy km s~ 


20(1) 


17(1) 


/cO(2— 1) 


Jy km s~ 


62(3) 


15(2) 


AVcoa-o) (FWHM) 


kms -1 


450(30) 


360(30) 


^CQ(2-1)' ^00(1-0) 




0.78(5) 


0.6(l) a 


^co 


kpc 


0.9(1) 


3.7(3) 


^major/^minor 




1.3(1) 


1.4(2) 


M gm 


M© 


4.2(2) x 10 9 


3.6(2) x 10 9 


M dyn 


M© 


4(1) x 10 10 


5(2) xlO 10 


M gas /M dyn 




0.11(3) 


0.07(3) 



CO luminosit i es ha ve been computed using equation 3 in 
ISolomon et al] (1997). Numbers in parentheses indicate estimated 
measurement errors in units of the last significant figures. a 
Corresponding to the central region (Rco = 1.1 kpc) of the SW galaxy. 



two lines. CO(1-0) and CO(2-l) spectra are shown in Fig. [2] 
and the position-velocity diagrams along the major axis of each 
galaxy are in Fig. [3] 

The CO distribution of the NE galaxy is highly concentrated 
(Figs.QJ andQJ). Given its compactness, it is not surprising that 
~100% of the CO(1-0) emission of this galaxy is recovered by 
the interferometer (Fig. |2^). Very likely most of the CO(2-l) 
emission is also recovered, as the £crj(2-i)/^co(i-0) luminosity 
ratio is ~0.78, similar to the typical lum inosity ratios found in 
normal spiral galaxies dBraine et alJll993h . 

The emission is resolved in both transitions, but it is in the 
CO(2-l) intensity map, due to the smaller synthesized beam, 
that the structure of the galaxy is more clearly seen. The centroid 
of the molecular gas distribution o f the NE galaxy co incides with 
its 15 GHz radio continuum core (Naga retal.120031) within 0.1". 
The deconvolved major axis of the CO emitting region (at zero 
intensity) is 1.7 kpc (1 .2"), s imilar to the sizes found in ULIRGs 
in advanced mergers (DS98). From the ratio between the major 
and minor axis of the CO distribution we derive an inclination 
for the galaxy of i ~40°. 

The symmetric CO(1-0) line profile and the resolved dis- 
tribution of the CO(2-l) intensity map (Figs. [2^ and QJ) hint 
that the molecular gas is distributed in a circumnuclear disk 
or ring. If such were the case, the asymmetric CO(2-l) line 
profile (Fig. [2}:) could be explained by a differential excitation 
or by some missing flux at the velocity range between to 
200km s -1 . The CO(2-l) position-velocity diagram along the 
mayor axis of the galaxy shows a linear increment of the ra- 
dial velocity as a function of the distance to the center up to 
Rco = 0.6" (Fig. [3^). The terminal rotational velocity at this 
distance from the nucleus is ~275kms _1 (cf. Figs. [2^ and [3^). 
Therefore, we derive a dynamical mass corrected from incli- 
nation of 4 10 10 Mo- The molecular gas mass in the same re- 
gion estimated from the CO(1-0) flux, and assuming a conver- 
sion f actor a - 0.8 M Q (Kkms -1 pc 2 )~' suitable for ULIRGs 
(IDS98I) . is M gas = 4.2 10 9 M , i.e. 1 1 % of the dynamical mass. 

The molecular gas distribution of the SW galaxy is more ex- 
tended than that of the NE galaxy, with a compact peak that re- 
mains detected even at the highest resolution. The deconvolved 
size of the major axis of the large CO(1-0) emitting region is 
7.4 kpc. From the ratio between the major and minor axis of 
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Fig. 3. a) CO(2-l) position-velocity diagram of the NE galaxy along the direction of its major axis (see Fig. QJl) at a velocity 
resolution of 30kms _1 . Contours for the CO intensity are 3cr to 13<x in steps of 2cr (lcr = 4.25 mJy beam -1 ). The black cross marks 
the central position and velocity of the NE galaxy as measured in the CO(1-0) map. b) CO(1-0) position-velocity diagram of the 
SW galaxy along the direction of its major axis (see Fig. QJ)) at a velocity resolution of 30kms _1 . Contours for the CO intensity 
are 3<x to 9<x in steps of 2<x (lcr = 1.57 mJy beam -1 ). The black cross marks the central position and velocity of the SW galaxy as 
measured in the CO(1-0) map. 



the CO distribution we derive an inclination for the galaxy of 
i ~ 45°. In Fig.[2j3 we compare the deduced single beam CO(l- 
0) spectrum of the SW galaxy with the CO(1-0) line emission 
detected by the PdBI. We estimate that -100% of the CO(l- 
0) emission of the galaxy is recovered by the interferometer. 
However, probably most of the extended CO(2-l) line emission 
of the galaxy is lost due to the higher frequency of the CO(2-l) 
transition. In fact the ^co(2-i)/-^co(i-0) l ummos ity rat i° calcu- 
lated in the central 2.2 kpc of the galaxy is ~0.6, a low value 
when compared to those found in spiral galaxies. 

The double horned CO( 1-0) and CO(2-l ) line profiles of the 
SW galaxy (Figs.|2]3 and[2jl) hint that the molecular gas is dis- 
tributed in a circumnuclear disk or ring. The position-velocity 
diagram along the mayor axis of the galaxy in the CO(1-0) tran- 
sition is plotted in Fig. |3p. The radial velocity increases as a 
function of the distance to the center up to Rco = 0.8" and then 
remains constant up to Rqo = 2.6" (3.7 kpc). The terminal rota- 
tional velocity is ~170kms _1 (cf. Figs.|2]3 and|3j)), therefore we 
derive a dynamical mass of 5 1O 1() M0 within this radius. The 
molecular gas mass in the same region is M gas = 3.6 10 9 Mq 
and M m /M dyn -7%. 



4. Discussion and conclusions 

From our CO(1-0) interferometric observations of IRAS 
10190+1322 we conclude that the two galaxies of the system 
have similar CO(1-0) luminosities, i.e. similar gas reservoirs. 
The main difference between them is the more compact molec- 
ular gas distribution of the NE galaxy, which is a consequence 
of the initial configuration of the encounter (see Murphv et al. 
1200 ll for a detailed discussion), being this galaxy the most 
affected by the interaction. The NE galaxy is cl assified as a 
LINE R and the SW galaxy as a HII galaxy by iRupke et all 
(2002 ]). Therefor e cont rary to the iDinh-V-Trung et al.1 d200Th 
and Evan s et ail (|2002) conclusions, we find that the active 



galaxy does not have the larger molecular gas amount, but a 
hig her concentration o f gas in its nuclear region. 

lFarrahet~ai]d2003l) obtained from starburst+AGN model fit- 
ting of IRAS 10190+1322 SED a total star formation rate of 
195 Mo yr~'. At this rate the total molecular gas reservoir of 
the NE galaxy will be exhausted in ~20Myr, considering that 
the bulk of the star formation takes place in this galaxy (see 
Sect. Q] and the discussion of our 2.8 mm continuum detection 
in Sect. This period is shorter than the time needed for the 
two galaxies to coalesce or experience a second close passage 
(> 50 Myr), assuming a transversal velocity of the order of the 
measured relative radial velocity (~140kms _I , see Table[TJ. The 
second encounter or final merger of the system may result in the 
concentration of the molecular gas in the central region of the 
SW galaxy and the restarting of the ULIRG phase with a new 
starbursting event. 
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